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Abstract 

A new type of perturbative expansion is built in order to give a rigorous derivation 
and to clarify the range of validity of some commonly used model equations. This 
model describes the evolution of the modulation of two short and localized pulses, 
fundamental and second harmonic, propagating together in a bulk uniaxial crystal with 
non- vanishing second order susceptibility x^ 2 \ an d interacting through the nonlinear 
effect known as "cascading" in nonlinear optics. The perturbative method mixes 
a multi-scale expansion with a power series expansion of the susceptibility and must 
be carefully adapted to the physical situation. It allows the determination of the 
physical conditions under which the model is valid: the order of magnitude of the 
walk-off, phase-mismatch, and anisotropy must have determined values. 

1 Introduction 

1.1 A methodological problem 

Many scientists, in recent times, came to the idea that physics is intrinsically nonlinear and 
that most of the linear theories have been introduced because they were mathematically 
tractable, contrarily to the nonlinear problems. Nowadays, the theoreticians of nonlinear 
physics almost always try to keep the largest part of the ideas involved by the linear theories 
of the physical problems considered. Both mathematical tools and physical concepts are 
conserved. The nonlinear theory appears then mainly as a perturbative approach in which 
the linear theory is the zero order. As an essential example, the Fourier theory and 
the plane monochromatic wave concept are the main tools in the linear theory of wave 
propagation. The theory of nonlinear wave propagation considers some wave packet, that 
is in fact a monochromatic plane wave (i.e. the linear wave) modulated by some slowly 
varying envelope (which is thus some first correction to the linear theory). This theory 
involves thus at least two space scales: that of the wavelength and that of the pulse 
length, which are treated separately by the mathematical model. The natural frame for 
the derivation of the nonlinear models is thus the so-called multiscale expansion method. 
From the mathematical viewpoint, formal computations in this frame can be performed 
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in a perfectly rigorous way. Convergence proofs have also been given ||. They enhance 
and clarify the validity of the formal computations. 

A multiscale expansion, as well as any asymptotic expansion, has a structural drawback. 
It involves some small parameter e, that is infinitely small in the mathematical theory, 
while it takes a small but finite value in the physical applications. Consider then some 
physical quantity, small in relation to the physical value of e, but independent of it. 
Because it is finite while e is not, the mathematical theory considers it de facto as if it 
were large in relation to e. This may sometimes introduce erroneous results, in so far that 
some further approximation is not taken into account. To introduce, beside the slowly 
varying envelope approximation and mixed to it, a further approximation of another type 
yields great technical difficulties. Such an approach mixing a multiscale expansion to an 
approximation of another type has already been done for the treatment of damping for 
soliton propagation, in particular for electromagnetic waves in ferromagnets [15, The 
damping parameter is identified to some power of the perturbative parameter e, and all 
involved quantities are expanded in a power series of the damping constant. 

In other situations, the ansatz must be constructed through a careful analysis of the 
physical situation. We study in this paper a situation in which such a mixed expansion is 
needed. It is presented from a methodological point of view, and should be of interest for 
the readers who intend to apply the same kind of perturbative approach to other physical 
situations. The proper physical consequences of the result will be published elsewhere [Fy 



1.2 Solitons through cascading 

The physical situation we consider is the propagation of optical solitons using the so- 
called cascading phenomenon, viz. second-harmonic generation and back-conversion to 
the fundamental, that arises in a dielectric medium with a non-vanishing second-order 
nonlinear susceptibility X ■ Several theoretical results have been published in order to 
show that this phenomenon could be responsible for the stabilization of a short pulse 



and for solitonic behavior [jy, |9L |1J, |lq| . These theoretical works study by analytic and 
numerical methods some partial differential system, derived in the frame of the slowly 
varying envelope approximation, and analogous in its form to the nonlinear Schrodinger 
(NLS) equation, at least for the linear part. The nonlinear part of this system gives account 
for the interaction between the fundamental and the second harmonic. Some applications 
of this model have been found in good agreement with experimental results. However, the 
way it is derived in the literature is quite phenomenological. It appears thus to be worth 
giving a rigorous formal derivation of that model, and to clarify its range of validity. 

Two waves propagating together are considered here, and their interaction is expected 
to stabilize the pulse. Therefore two conditions must be satisfied, at least in an approxi- 
mate way. First the interaction must be resonant, thus the waves must be phase-matched. 
If not the interaction will not be strong enough to yield the expected effect. Second the 
two waves must propagate together during a long enough time. Thus their group veloci- 
ties must be equal, or at least very close together. We determine precisely the order of 
magnitude of the allowed small deviation from phase-matching and small walk-off. We 
restrict for technical reasons our attention to crystals having some particular symmetry. 
However, the required phase and group velocity matching conditions cannot be realized in 
any crystal of this type: the linear dispersion relation of the material must satisfy several 
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conditions. The order of magnitude of some physical quantities must be comparable to 
powers of the perturbative parameter characteristic of the multiscale expansion. 



2 Analysis of the physical problem 

2.1 The equations for a uniaxial optical crystal 

We consider a bulk crystal with a non-zero second order susceptibility (x ) tensor, opti- 
cally uniaxial from the linear point of view and perfectly transparent. These hypotheses 
are satisfied by potassium dihydrogen phosphate (KDP) and analogous materials or by 
lithium niobate, as examples. Wave propagation in such a medium is described by the 
Maxwell equations that reduce to 



V A (VAE) = -^3} (e + 3) , (2.1) 



where c is the speed of light in vacuum, and 

D = x (1) * E + x {2) * EE 
ft r t' 

x (i) (t _ t')E{t')dt' + / x (2) (t-t',t- t") : E(t')E(t")dt"dt'. 

) J —oo J —oo 



(2.2) 




and x*- 2 -* are respectively the linear and the quadratic susceptibility tensors Q. We 
denote by x^( u ) the Fourier transform of X^(^) so that x^W = /_ °° X^i^e^dio. 
Then, for a uniaxial crystal, the coordinate frame being chosen so that the optical axis is 
the z-axis, the linear susceptibility tensor satisfies 

77^ 

X (1) = ( ^ I . (2.3) 

n (uj) and n e {uj) are called respectively the ordinary and extraordinary indices [Q]. A wave 
propagates along the optical axis with the velocity c/n Q . If the propagation direction makes 
some nonzero angle 8 with this axis, two waves can propagate, with different speeds and 
polarizations. They are referred to as ordinary and extraordinary waves. The charac- 
teristic feature of the extraordinary wave is that its group velocity is not parallel to its 
phase velocity, which never happens in isotropic media. If 9 = vr/2, the wave polarized 
perpendicular to the optical axis is the ordinary one, the wave polarized along the axis is 
extraordinary, and propagates with the velocity c/n e . 



2.2 The phase-matching angle 

No wave interaction can occur far from the phase resonance, or phase-matching. Therefore 
we must work close to it. We denote by v v (uS) and v ip (2uj) the phase velocities of the 
fundamental and second harmonic respectively. The phase-matching condition obviously is 



v v {(jS) = Vtp{2tjj). 



(2.4) 
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It is well-known [17] that this phase-matching condition can be realized by some adequate 
choice of the value of the angle 6 between the optical axis (z-axis) and the propagation 
direction. Several situations arise, depending on the relative magnitude of the indices in 
the considered material and on the polarization of both waves. We consider only one of 
these cases, of major importance for applications. For KDP and many other materials, 
the inequalities n {2uj) > n {uj) and n (uj) > n e {uj) are satisfied. It is easily shown that in 
this case the phase-matching angle exists when the fundamental wave is ordinary and the 
second-harmonic extraordinary, but does not exist for other polarizations, except if both 
waves are extraordinary. In the latter situation, the experimental values of n e and n Q do not 
yield a real phase-matching angle value. We restrict thus the study to the particular case 
where the fundamental wave is ordinary and the second-harmonic extraordinary, which 
is the only physical situation in which the phase-matching angle exists for the materials 



under consideration. The phase- matching condition is then [17] 



sin 



n- 2 H-n- 2 (2oQ 
rie 2 (2uj) -no 2 (2w)' 



(2.5) 



The first condition to be satisfied by our ansatz is that the angle 6 between the propagation 
direction and the z-axis must satisfy the phase-matching condition ( |2.5[) , at least in an 
approximate way. 



2.3 Transverse velocity 

The dispersion relation uj{k) follows from fl2,3| ), the value of the group velocity is then 
computed according to v = doo/dk. While the group velocity of the ordinary wave is par- 
allel to the phase velocity, that of the extraordinary wave has a transverse component [17]. 
It is 

c 2 k ( 1 1 



with 

.2 , 2,2 



cos 2 6 n' sin 2 6 n[ 



A = or + c l k l I "° + : " e (2.7) 

V n i n l ) 

(k is the wave vector and u the pulsation of the fundamental, n' Q = dn /dtu). Because 
we expect that the walk-off is small, v ex must be close to zero. This can be achieved in 
two ways: either by setting n e close to n , or when sin 26 is small. Try first to satisfy 
the former condition: anisotropy must be weak. However, anisotropy is needed for the 
phase-matching that must be realized at least in an approximate way. The difference 
(n e — n ) appears indeed in the denominator of fl2.5|) , and therefore it cannot be zero. The 
anisotropy, or the difference (n e — n D ), must be small but not too small. This is the second 
condition imposed to our ansatz. It can be satisfied in many real materials. 

It is likely also necessary that sin 26 is small. We note in fl2.5|) that, when n e (2uj) = 
n (tu), then sin 2 6 = 1. In this case 6 = ir/2 and sin 26 = 0. We deduce this way 
a third condition to be satisfied by the ansatz: n e (2uj) must be close to tl (uj). Thus the 
achievement of phase-matching and reduction of walk-off necessitate that n (2uj) and n e (u) 
are very close together, and that anisotropy is weak. The ansatz used below is built in 
order to satisfy these conditions. It is seen from published experimental data Q that they 
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are realized in many real materials which are in fact the most commonly used for practical 
applications of the second harmonic generation. Further it was very early recognized that 
the arising of the transverse group velocity or walk-off is the main limitation to an efficient 
second harmonic generation at phase- matching angle ||, and the fact that the condition 
"n e (2u;) close to n (o>)" ensures the reduction of the walk-off in this situation follows 
straightforwardly from the properties of light propagation in a birefringent medium. 



3 An adequate ansatz 

3.1 The scaling parameter 

As usual in the frame of multiscale analysis, we introduce some small parameter e. It 
characterizes the pulse size in the following way: we call L the common order of magnitude 
of the pulse length and width, and A the wavelength. Then the perturbative parameter is 
defined as 

(3-D 

The electric field E is expanded both in a power series of some small parameter e, and in 
a Fourier series of some fundamental phase <j) as 

E= £ l Efe ip4 '. (3.2) 

l>l ,V& 

Iq fixes the order of magnitude of the main term. It is discussed below. The amplitudes Ef 
are functions of slow variables that are defined below. The main slow variables characterize 
the pulse shape and its propagation at the group velocity. According to ( |3.1[) , they are 
variables of order e, defined by 

£ = ex, rj = ey, ( = ez, r = et. (3-3) 

We seek for the evolution of the pulse shape through propagation distances long in relation 
to the pulse length L, say L/e, as in the standard nonlinear Schrodinger model. On these 
grounds and because the propagation direction is close to the x-axis (see below), the slow 
variable giving account for the evolution can be chosen as 

£ 2 = e 2 x. (3.4) 

3.2 The carrier 

The phase <j) determines the carrier monochromatic plane wave. It thus determines the 
propagation direction. Recall that, because phase-matching must be at least approxi- 
mately realized, the angle 9 between the z-axis and the propagation direction must be 
close to the value given in ( |2.5D . On the other hand, it has been shown in subsection 2.3 
that the latter value must be close to tt/2. We write it as 

e = l~e a l, (3.5) 
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where 7 is a free parameter at this point, and the exponent a is positive and is determined 
below. The phase (j> is 

(p = k ■ x — cot, (3-6) 



with k = /c(sin#,0,cos#). Inserting the definition ( |3,5[ ) of 9 into the expression ( |3.6D of 
we would get 

<j> = kx + A;7£ a 2 - wt + O (e 2a ) . (3.7) 

A dependency with regard to some slow variable e a z appears (a > 0). For a > 1, this 
dependency can be incorporated in the amplitude, that depends on £ = ez. Thus only 
values of a less than 1 are interesting. On these grounds, and taking into account the fact 
that the perturbative expansion will involve only integer powers of e a , it is seen that the 
correct choice for a is a = 1/2. Further, due to the anisotropy, the wave vector norm k 
depends on the angle 6. Thus, according to (3.5), k depends on e. As e tends to 0, k must 



be expanded in a power series of \fe. A term involving \fex appears then in (^7|) . The 
terms of order e 1 and higher are incorporated in the amplitude. The complete phase (j> 
can thus be written as 

4> = k (l + a\fe) x + \fek^z — tot. (3.8) 

In the same way as k, all quantities such as group velocities or polarization vectors, that 
depend on the angle 6, are expanded in a power series of \fe. Therefore we must introduce 
an additional slow variable of half-integer order, 

ix = e^x. (3.9) 

The complete set of slow variables on which the amplitudes Ef depend is thus 
(£, r), C, r, £i, £2). They are defined by Q, Q, and Q. 

Leaving aside the half integer powers of e, this is mainly the "classical" multiscale 
expansion that leads ordinary to the nonlinear Schrodinger-type (NLS-type) models ([0], 
p. 495 sq.) describing long-distance propagation along the x-axis. Recall indeed that 
it is not necessary to introduce a priori the statement that the field is at first order a 
function of (x — Vt). Else, the linear transport equation, that is the equation giving 
account for the modulation propagation at group velocity, is obtained as a solvability 
condition of the perturbative scheme. The main originality of the expansion ( |3.2j )-( |3ljD is 
the introduction of an intermediate scale X\Je, z^/e in the expression of the phase <j). The 
expansion therefore involves not only 3 scales for the longitudinal space variable, as in 
the usual NLS-type expansions, but 5 different scales (including £1). The term y/ek^z in 
expression ([T^) for <p represents a deviation for an angle yfe^ from the propagation along 
the x-axis. 



3.3 Expansion of the susceptibility 

Here comes the difficult point which is the introduction of a second approximation beside 
the multiscale expansion. It has been seen that the components of the linear susceptibility 
tensor must satisfy several hypotheses that allow certain approximations to be made. The 
problem is that the only consistent way of describing mathematically these hypotheses 
requires the tensor to be written as a function of the perturbative parameter e. 
However, e represents the pulse size, and does not depend on it from the physical 
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point of view. Nevertheless, we assume that the linear susceptibility tensor x^ can t> e 
formally expanded in a power series of y/e in the following way: 



1 + x « = n 2 + y/iffi + ex? + ^X^ + • • • • (3.10) 

2 2 



n 2 is an isotropic index, and the fact that lim (l + X ) = n2 gives a partial account for 
the second condition of § 2.3, the weak anisotropy hypothesis. A particular definition of 
the coefficients Xj can De given. It is of no interest for our purpose, and is left for further 



publication |13 |. 



3.4 The amplitude 

An important difference between the ansatz ( |3.2|) and the standard one for NLS is the 



order of magnitude of the leading term, that is the value of Iq in (3.2). It is determined 
as follows. The model that we intend to derive and to justify is used in |l6|, 14, ||. It has 
the following form: 

idfr<p + Dl<p = A 1 ipip*, (3.11a) 
id^ + D ,2 ^ = A 2 ip 2 , (3.11b) 

where Di and D'^are second-order partial differential operators relative to the slow variab- 
les (r — £/V), r], ( that describe the shape of the pulse. A\, A 2 are some constants. 
cp and tjj are the amplitudes of the fundamental and second harmonic respectively. These 
amplitudes are assumed to have an order of magnitude e °. Furthermore, £, 77, C> r have an 
order of magnitude e, and ^2 an order e 2 . Therefore the left-hand-side of equation ( [3.11 ) 



has the order e 2+l ° , and its right-hand-side has the order e 2l ° . We must thus have 2 + Iq 
21q, and Iq = 2. 

This differs from the case of the cubic NLS equation, which is 

idzJ + Dlf + Af\f\ 2 = 0, (3.12) 

with analogous notations, / being the modulated amplitude. There we would have Iq + 2 = 
3/o, thus Iq = 1. This has a very important physical meaning: the present model requires 
input energy pulses from a smaller range than the NLS solitons obtained either in cubic 
media by the Kerr effect or in quadratic media by cascading, far from phase-matching || . 
This is due to the fact that the interaction is resonant, or very close to resonance. Note 
that the multiscale analysis justifies this difference between the required power inputs 
without making any special assumption about the values of the nonlinear susceptibilities. 

4 Solving the perturbative scheme 
4.1 A technical point 

The order by order resolution of the perturbative scheme can now be performed. The 
detail of this derivation is given in the appendix. The vector amplitudes whose evolution 
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is under investigation are polarized so that 

E\ = [ <p V E\ = I I 11.1) 





The transport equation giving the group velocity is expected to appear at order e' 0+1 = e 3 . 
We find that the amplitude ip of the fundamental satisfies the partial differential equation 

df(p + -d T (p = i5\(p. (4.2) 

v 

The constants 5\ and v are specified below (equation (p~g|)). The solution of equation ( [1.2| ) 
is 

V&Zu&T, V X) = w Ui,&,T-^,ri,() e i6 ^. (4.3) 

It describes both a propagation at speed v and a linear phase modulation measured by 
the coefficient Si. 

The expressions of the constants v and Si are given by solving the multiscale expansion. 
They can be interpreted, and written in a convenient way, by using some power series 
expansion of the wave numbers. We introduce the following notations: 

h 2 s = n 2 + ^Anl + ex ( ±l + e^x ( p -\ for s = o or e, (4.4) 

' 2 ' S 

k a = k Q (uj) = LOn °^ = j~ + + ski )0 H , (4.5) 

k e = k e {2oj) = 2uJUe ^ = koe + y/iki + £ki }C + ■■■ . (4.6) 

C 2 



For the derivatives of these quantities, the notations used are 

±k p , {u) and k' p>e = ^ 



k'p,o = -l-kpA 10 ) and k p,e = Ail, A k P,ei 2uJ )- ( 4 - 7 ) 



Then the values of the coefficients are found, 

<5i = fci, - fc °'° 7 and i = A;q i0 . (4.8) 

Thus the phase factor e lSl ^ accounts for the corrections to the x-component of the wave 
vector k at order e, in norm (term ki :Q ) and in direction (2 nd term). Similar corrections 
appear at following order e 3+ 2, it is the factor e lp1 ^ 1 in formula (|A.9|) , and for the second 
harmonic: the factors e lp1 ^ and e %P2 ^ 1 in formulas ( 4.13a ) and (|4.13b[) . These corrections 
are a characteristic feature of the kind of "mixed" expansion used here. 

Note that it is decisive, for going further in the perturbative expansion, that equa- 
tion ( |4.2| ) has its particular form. The multiscale expansion frame considered describes 
indeed the long-distance evolution of some state, that is a steady state in the first ap- 
proximation. Usually, the order zero is a plane wave, invariant in amplitude, thus order 
one can be considered. The order one presents Galilean invariance, thus order two can 
be considered. Higher order make sense mainly if order two describes a soliton, which 
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is another kind of invariance. Therefore, the arising of a ^-dependency should a priori 
forbid to go to the next order, except if this dependency presents some invariance. As an 
example, an exponentially decreasing ^-dependency, describing absorption, would reduce 
all following orders to small linear corrections to this decreasing, that would give account 
for the main physical behavior. There is, however, some invariance here, because the 
^-dependent factors are unitary phase factors, and we can pursue the expansion. 

4.2 The phase-matching condition 

The first order of the dispersion relation, found at first nonzero order e 2 , is 

c c 

Because k is unique, the first order of the phase-matching condition is obtained as 

n(1uj) = n{u). (4.10) 

The common value of n(u) and n(2u) is written n below. 

At order e 2+ 2, the anisotropy appears through a condition involving ak, the term of 
order yfe in the expansion of k. We choose, as written above, the fundamental and the 
second harmonic as ordinary and extraordinary waves respectively, a takes the value 

a = ^M = MM. (4.11) 

2n 2 2n 2 v ; 

Because a is unique, the second order of the phase-matching condition is got at this point. 
It is 

An 2 (w) = An 2 (2w). (4.12) 

Thus the two first orders of the phase- matching condition are a direct consequence of the 
ansatz, precisely of the unicity of the wave vector in the phase 4> defined by ( |3.8| ). A more 
complicated ansatz involving two fundamental phases could be envisaged. However, it 
would not be able to give account for the wave interaction. It is indeed proved below that 
the following terms in the expansion of the linear susceptibility must also coincide 
(see equation ( |4.16| )). The same kind of considerations would prove the same result for 
first order, if it was not set a priori in the ansatz. 

The following order of the phase-matching condition appears at order e 4 when wri- 
ting down the nonlinear evolution equation for the complex amplitudes ip and ip of the 
fundamental and second harmonic respectively. They are equations ( 3.1 1| ) , while the 
amplitudes ip and ip satisfy 

i p = i p 1 (^,r-^-^, V ,C-^i) e^+P^\ (4.13a) 

\ V V 1 J 

= Vi t - £ - ^, v, C - rti) e'M+nM. (4.13b) 

Equations ( |3.11| ) are expressed in a coordinate frame that moves with the wave, at a speed 
equal to the group velocity of the wave up to order yfe. The coordinates are defined by 

t = t , fj = T), C = C-7£b 6 = 6- ( 4 - 14 ) 

V Vl 



Mixed Perturbative Expansion: the Validity of a Model for the Cascading 



135 



In this frame, there is no dependency of the amplitudes relative to the variables £ and £1. 
On the other hand, phase factors depending on these variables appear. These factors are 

exp i [(62 - 25i)£ + ( P 2 - 2pi)6] (4.15) 

in the evolution equation for cp, and the inverse in the evolution equation for ip. These 



two equations yield equations ( 4.21a ) and ( 4.21b ) respectively, after cancellation of the 



factors mentioned. This yields two approximate phase-matching conditions, in addition 
to ( 4,11| ): 62 = 26i and P2 = 2p\. The first condition can be written as ki >e = 2fci )0 , or as 



X W(2a;)=xi 1) M. (4-16) 



4.3 The phase-matching angle 



Conditions ( |4.1C| ), ( |4.12 ) and (4.16) do not involve a phase- matching angle. They state 



that the approximate phase-matching must be realized "spontaneously" by the material up 
to order e. The last of the phase-matching conditions obtained in the previous subsection is 

P2 = 2pi. (4.17) 
It can be satisfied by some particular choice of the angle y/e'y according to 

„„ ks J2u) - 2ks (uj) 

7 2 = — 2'° . (4.18) 

7 lo An 2 (u;) - An 2 (2w) 1 ' 

Formula ( |4.18| ) coincides with the usual phase-matching condition ( |2.5| ), taking into ac- 
count the power series expansion of x an d the "spontaneous" approximate phase- 
matching, up to order e. A deviation from the exact phase-matching of order e 2 is possible. 
It is described by the terms Bnp in equation (4.21a) and B21P in equation ( 4.21b| ). The 



present multiscale expansion proves that the model is by no means valid for a larger de- 
viation. In such a case, only a model involving non-resonant interaction can describe the 
physical phenomenon. Completely integrable models such as the nonlinear Schrodinger 
equation (NLS) or the Davey Stewartson system (DS I and II) can be derived in this frame, 
but involve much higher intensities: of order e instead of e 2 in the present case [|ll], 12]. 



4.4 The group velocity matching 

The phase-matching is a phase velocity matching, but for short pulses a group veloci- 
ty matching is also necessary. Without it, the two pulses cross each other without ha- 
ving time enough to interact. Mathematically, the group velocity matching condition 
appears through the requirement that the variables (f, fj, C, £2) are uniquely defined by 
equation ( |4. 14 ) . In other words, the frame moving with the pulses, in which equation ( p. 11 ) 



describes the evolution of the waves, must be defined in a unique way, thus must be the 
same for both waves. 

At order e 3 , the first order term of the group velocity is found. It is 1/v = k' Q for the 
fundamental and 1/v = k' 0e for the second harmonic. According to the above statements, 
the velocity v is the same for both waves. If these velocities differ, it is found at next 
order that no interaction between the waves is possible. From the formal point of view, 
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the £1 and ^-dependencies of tp and ip would be linear if both group velocities have not 
the same value. For this reason, we must have k' Q = k' Q : this means that the group 



velocities coincide at order zero. Using the zero order phase-matching condition ( [4.10 ), 
this condition can be written in terms of the index n as 



n'{io) = 2n'(2w). 



(4.19) 



The first correction l/v\ to the inverse of the group velocity is obtained at following 
As at previous order, this correction is l/v\ = k'± for the fundamental and 

9 :° 



order e 3+ 2 



1/vi = k\ for the second harmonic. Therefore k\ = k\ or, in terms of the indices, 



2(A^)'(2a;). 



(4.20) 



An explicit walk-off appears in the final equations ( f4.21a| )-( |4.21b| ) through the terms in- 
volving a first order r-derivative. It is measured by the difference between k[ and k[ . As 
the deviation from phase-matching, the walk-off has a fixed maximal order of magnitude, 
that of e. If this group velocity matching is not satisfied, the model yielded by equations 
(l4.21a| )-( [Olb|) is theoretically no more valid. This means physically that for a pulse size, 
a power input, and a propagation distance in accordance with the scaling, the two waves 
do not coexist at the same place during a long enough time, and cannot interact. 



4.5 The asymptotic model 



As mentioned above, the evolution equations of the form (3.11) for the complex amplitu- 
des (p and ip of the fundamental and second harmonic respectively are obtained at order e 4 . 
They reduce to 



2i/e% - B l + 2ik ( k' ho + ^ ) d T 



2ikj 



- kk"d 2 T + dl + d}- 2ikcryd c 
Aikd^ -B 2 + Aik (k' l e + Q d T 
2An 2 + An 2 e ) <9 C - 2kk^ e d 2 T + d 2 + 9% 



if> = A 2 ip 2 



(4.21a) 



(4.21b) 



(The notations An 2 , kj jS (s = o, e) are defined by equations (|4.4j) to (^4.6[) ). The coeffi- 
cients, 2ik or 4ik of <9g 2 , and — kk" (= — fco,o&o o) or — 2fcfco e of d 2 , are the usual coefficients 
in the NLS models in nonlinear optics. A transverse velocity term, proportional to d^cp, 
and a longitudinal one, proportional to d T cp, appear (the coefficient a is given by equa- 
tion ( [4, ll| )). k[ and k[ are the corrections to the inverse of the group velocities at 
order e. 

The coefficients B\ and B 2 are some corrections to the phase vector at this order. They 
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have the expressions 
Bi = 



-2kk 2 ,o + 



3(An^) 2 ( 



Bo 



4n 4 
Akk 2 e + 4k 2 



k 2 

7 2 + T 7 4 , 



(4.22) 



L((An2) 2 + 4A^An e 2 -8(A^) 2 



8n 4 V v 



7 2 + fe 2 7 4 - 



(4.23) 



The nonlinear interaction constant A\ and ^2 are defined in the general case by equations 
( A.24| )-( [A.30| ) (in appendix). Using the symmetry properties of the x^"tensor, inclu- 
ding the so-called complete symmetry property, we have A 2 = 2A±, at least for the 42m 
and 3m symmetry classes, for which the common value of these coefficients is given by 
equations (|A.27|) and (|A.29| ). 

Equations ( 4.21a ) and ( 4.21b| ) yield the sought model. They are analogous to the 
equations given by Menyuk in []l4| , but here in a complete (3 + 1) -dimensional version. 
Kanashov and Rubenchik shown || that no collapse occurs in this model, and that a soliton 
solution exists in the sense of a stable localized pulse. System ( 4.21a| )- (|4.21b ) is written 
in a coordinate frame that moves with the wave at a speed equal to its group velocity up 
to order yfe. It is seen that the group velocity of both fundamental and second harmonic 
must coincide up to this order. Walk-off terms of order e appear in equations (4.21a) 



and ( 4.21b ). These are the terms involving first order r- and (^-derivatives. In an analogous 
way, it is seen that the approximate phase-matching must be realized "spontaneously" by 
the material up to order e, and once more at order e 3//2 by some particular choice of the 
angle y/e^y, according to equation ( f4.18| ). 

Consider now the expansion Q3.10| ) of the susceptibility tensor, that also has the ex- 
pression ( |2.3| ). Because we must have an isotropic zero order term, \fe is the order of 
magnitude of the anisotropy. The interaction must be phase-matched "spontaneously" up 
to order e, thus the order of magnitude of the difference (n (u>) — n e (2iv)) must be e 3//2 . 
Last, the group velocity of both waves must be identical up to order yfe. Thus their 
difference must have an order of magnitude e. Finally, the matching conditions are 



nl(uj) 



eo(Vi) 



G 0(e), 



n 2 (w) - 
dn 2 Q (uj) 

du> 
n 2 M - 



is the fundamental frequency and e the perturbation parameter. 



dn 2 e (2uj) 
dio 



n 2 (2w) G O (e 



J/: 



(4.24a) 
(4.24b) 
(4.24c) 



5 Conclusion 



A long-distance propagation model (equations ( 4,21a| )- Q4.21bp ) for two interacting waves 



with different velocities has been derived through a rigorous multiscale analysis. It uses 
a new type of expansion that mixes a multiscale expansion, and an expansion in a power 
series of some physical quantity, that here represents mainly the anisotropy. This expan- 
sion involves 5 longitudinal space scales, instead of 3 in the usual NLS-type models. The 
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supplementary scale is related to approximate conditions satisfied by the linear suscep- 
tibilities. The expansion considered will be useful for the mathematical justification of 
asymptotic models of the above mentioned type. This rigorous derivation has allowed the 
determination of the validity conditions of the model. 

From a more concrete point of view, the model equations (4.21a-4.21b) describe the 
evolution of the modulation of two short localized pulses, fundamental and second har- 
monic, propagating together in a bulk uniaxial crystal with non-vanishing second order 
susceptibility x^ 2 \ an d interacting. Much theoretical work has already been done about 
this system. It has been proved that no collapse occurs ||, analytical solutions have 
been found in (1 + 1) dimensions [pf l, and numerical simulations have shown the existence 
of stable solutions [16|. We have determined the physical conditions under which the 
model is valid: both phase and group velocity of the waves must be close enough together. 
Group velocities are matched when anisotropy is small. Thus the phase-matching angle 
must be small. Therefore, the phase-matching must be realized by the intrinsic properties 
of the material and for any propagation direction with some precision. Hence the model 
cannot be valid in any material and at any frequency. X-Mozilla-Status: 0000 



Appendix. Details of the derivation 
A.l Dispersion relation 



In this appendix, we give the rigorous derivation of system ( 4.21a )-( 4.21b ), using the 



perturbative scheme described in Section 2. The convolutions involved by expression (2.2) 



of D are treated using the formulas derived in 11]. We solve it order by order: the first 
nonzero order, e 2 , gives that E^ s is zero for all s = x,y,z, and all integer p, except for 
p = ±1, ±2, and s = y,z, if the dispersion relation (|4.9| ) is satisfied. 

At order e 2+ 2 ; a first correction to the dispersion relation is obtained through the ex- 



pression ( 4.1 1| ) of the coefficient a. The polarization vectors E\ and E\ are also found, they 



are defined as (|4.1|). 92 and ip are the amplitudes whose evolution is under investigation, 
and the other E^ are zero (except for the symmetric terms p = —1, —2). 

A. 2 Order 3 and group velocity 

At order e 3 , the group velocity appears. Precisely, we get as a solvability condition for ip 
the equation $L%), which admits the solution (|4.3|) . In an analogous way, the ^-dependency 
of the second harmonic is obtained: 

^ = ^o(^,^,t-^,7 1 ,c] e 1 ^, (A.l) 



V 



with 



<5 2 = A;i, e -^£2! and - = k' 0>o . (A.2) 



Other conditions give 

E\ +h = I / I , El +h =\ ) . (A.3) 
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The functions / and g have to be determined. 
A. 3 Order e 3+ ^ 

The following order e 3+ 2 gives, in a similar way, the term of order yfe in the group velocity. 
The following equation is then obtained: 

2ik 

—2ikd^f — 2ikad^ip — 2ikd^(p = 2k5\f + -^-d T f 

+ 4*?> + % ((An 2 )'o, + 2Anl) d T p. (A.4) 



Equation (A.4) is solved as follows: using the change of variables 

f = r-i e = c, (a.5) 

V 

noting that ip = (po(f)e lSl ^ (equation (fOl)), and setting 

f = fo(ir)e i5 ^ (A.6) 

in equation (|A.4| ), it is found that c^/o expresses in terms of (fo and of its derivatives. 
Thus, relative to £, d^f is some constant A, and 

fo = Ai + B, (A.7) 

B being some constant. Because fo must be bounded as £ tends to infinity, A must be 
zero. This implies that fo is a function of f alone, and this yields some partial differential 
equation to be satisfied by ifo (or <p). This equation is obtained from equation ( |A.4| ) by 
simplifying the terms containing /. It is 

1 

%<A) + 7<9c ( A) H d T if - ipupo = 0. (A.8) 

Vi 



The solution of equation ( [A.§| ) is 

¥>o = ¥>i (&,T-^-^, V ,(-rti) e ip1 ^. (A.9) 
V v n ) 

1/vi is the first correction to the inverse of the group velocity 

— = k\ . (A. 10) 

pi is the next correction to the wave vector k, 

Pi = fc§,o + ^ A »oH7 2 - (A.H) 

The dependency relative to (C — 7£i) does not represent a walk-off: it simply gives account 
for wave propagation in a direction making the small angle jy/e with the x-axis. 
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In an analogous way, the ^-evolution of the second harmonic is found, 

9 = 9o { r 



with 



and 



1 = ^ e > 

V\ 2' e 



p 2 = ks 



W7 
2nc 



(An^(2o;) - 2An 2 (2uj)) . 



(A.12) 

(A.13) 
(A.14) 



Equation ( 4.13b| ) is deduced from equation ( A.12|) . As previously, the variables must be 



the same in expression ( |A.9| ) for <p, as in (4.13b) for ip. In particular, the correction 
term \ jv\ to the inverse of the group velocity must be the same in both cases. Thus 



k\ 



■O 



k' le . 

2 ' e 



(A.15) 



Other conditions give 




Ei 



El 



-19 + + & [2An2(2a;) - An*(2w)] ^ 


G 



(A.16) 



(A.17) 



F, G are functions to be determined. The expressions of E '\ j E^' x 1 are also obtained, 

d+2 J+2 



they are use in the remainder of the computation: 



E 



hx 



l 9 vf ~ ^2 d v<P> 



E 



2,x 
^3+1 



2n 2 
_7_ 
2n 2 



(A.18) 



-fG + -d ( g + ^ [2An 2 - An 2 e ] g 



+ 



7^ 



Z7 



2fc 



4/cn 



(A.19) 



(All indices in equation ( A.1S| ) are taken at 2u>). 



A. 4 The evolution equations 

At order e 4 , the evolution equations for tp and ip are obtained as follows: the equation of 
order e 4 for the y-component of the fundamental is a partial differential equation involving 



F, f, and <p. Using the variables £ and f defined by equation ( A. 5 ), and writing F 



Eo(£, f)e lSl ^, it is found that dgFo is equal to some expression depending on tpo, fo, and 

their derivatives. This expression does not depend on £. Thus, as shown above for /o, .Fo 
is a function of f only, and the terms depending on F cancel in the equation. 
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Then considering equation (|A.9|) , /o is 



/o = h Ui, 6,t-^,17,C-7&1 e^ 1 . (A.20) 



The notation holds for the partial derivative relative to £1 in the variables defined by 
equation ( A,20| ). Transferring (A.9) and ( |A.20 ) into the equation, the following result is 



obtained: d^fi is equal to some expression that depends on ip\ only. Thus d^fi does not 
depend on £1. Because f\ must be bounded as £1 tends to infinity, f\ is independent of £1 
too. Then f\ disappears from the equation, that reduces to the evolution equation (4.21a) 



for ip, taking into account the nonlinear term computed in the following subsection. 

The equation for the second harmonic is obtained in a similar way. It comes from 
the equation of the multiscale expansion at order e 4 , for the z-component of the second 
harmonic. In the same way as for the functions F and / in the case of the fundamental, 
it is seen that the functions G and g must have the form 



G = G It - J J e id2 S (A.21) 
and 

9 = 91 U2, T ~ ~~ C lb) e'&Z+nV. (A.22) 
Then all terms involving either G or g cancel in the equation, and the evolution equation 



( 4.21b ) for ip is obtained. 



A. 5 Computation of the nonlinear terms 

The nonlinear term in equation (4.21a) for ip is 



uJ 



2 



Pl' y , (A.23) 



where P^ ,y is the y-component in the coefficient P4 of e 4 e l< ^, in the expansion of the 
nonlinear polarization '■ EE. It is straightforwardly seen that 

Pi = 2 X {2) (2u, -u) : E\E 2 \ (A.24) 

The value of this nonlinear term depends on the symmetry properties of the %( 2 )-tensor. 
For the 42m class of crystals, to which belongs e.g. KDP, in the crystal axes frame Q, the 
non- vanishing components of the -tensor are 

(2) = (2) (2) = (2) (2) = (2) , A 25) 

Axyz A-yxzi Axzy Ayzxi Azxy Azyx- yr^.^u/ 

and all other x^-components are zero. Thus, if the axes of the coordinates frame are the 
crystal axes, because E\ and E\ are parallel to the z- and to the y-axis respectively, P^' v 
would be zero. Thus the crystal axes must be rotated for some angle (5 around the optical 
axis (the z-axis). We call (x',y',z') this rotated crystal axes frame. Relative to it, the 
nonlinear term becomes 

-PI =A 1 { I VV\ (A.26) 



uJ 2 



C 2 



1 
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with 



2oj 



2 



Ai = ^ sin 2/3 X%A 2 "i -w). (A.27) 
For the 3m symmetry class, the nonzero components of the X(2)~tensor are j|] 

(2) = (2) (2) (2) (2) (2) (2) 

A,xzx /^yzyi A,xxz A.yyzi A.zxx Azyyi A.zzzi 

(2) = _ (2) = _ (2) = _ (2) (A 2g) 

A-yyy /^yxx A.xxy A,xyx' V 1 *" 1 -'/ 

The coefficient A\ is independent from the angle introduced for the 42m class, and has 
the expression 



2uj 2 



M = r X£4,(w,o/). (A.29) 

Other values for this constant are obtained with other crystal symmetries. They are easily 
computed, or can be found e.g. in ||. 

In equation ( [4.21b| ) for ip, the nonlinear interaction constant Ai is defined in the general 
case by 

_^ P V = A 2 ^ 2 . (A.30) 
For both the 42m and the 3m symmetry classes, the use of the %^ -structure (A. 25) 



or ( [A. 28 ), and the complete symmetry property of the x^"t ensor ) proves the equality 



A 2 = 2A X . 

This completes the derivation of the model equations ( 4.21a| )-( [4.21b ). 
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